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Tliis paper is based on the presentation made at the XIX th URSI General Assembly, Helsinki, 
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, A REVIEW OF RECENT RESULTS ON SPREAD F THEORY 


I. Introduction 

Spread F, as exhibited by diffuse echoes on ionograms, was 
discovered some forty years ago by Booker and Wells (1938). Until 
recently, there has been only much statistical data concerning 
Spread F. However, in the last few years with advances in radar 
backscatter measurements, in situ measurements, and theoretical 
and numerical simulation techniques a clearer picture of the funda¬ 
mental plasma instability mechanisms causing equatorial Spread F 
phenomena has been evolving. It should be emphasized that this paper 
will deal only with equatorial Spread F theory and, in addition, dwell 
only on those theories using plasma mechanisms as a basis. For 
experimental results, the works of Farley et al. (1970), Dyson et al. 
(1974), Kelley et al. (1976), Woodman and La Hot (1976), Morse et al. 
(1977), and McClure et al. (1977) should serve as good references for 
the interested reader. 

In speaking about Spread F one should at least show its basic 
manifestation on an ionogram (see Fig. 1). After all, the terminology 
Spread F emerged from the results of ionosonde traces such as those 
exhibited in Fig. 1. The multiple traces are caused by magnetic field 
aligned irregularities. If there were no irregularities, a single 
crisp trace rising to the right would be exhibited on the ionogram. 

At this point we want to remember the basic equatorial geometry. 

Figure 2 exhibits the basic equatorial nighttime ionospheric F region 


geometry, i.e., the geometry under which equatorial Spread F occurs. 
Note: Manuscript tubmitted November 16, 1978. 









N(y) represents the background electron density as a function of 
altitude (y). Gravity, £, points down, the ambient magnetic field, B, 
is horizontal (pointing north) and k represents a horizontal perturba¬ 
tion (in the westward direction). The maximum in the electron density 
profile is the F peak. The underside of the profile steepens at night 
due to chemical recombination effects and electrodynamic forces. The 
E region has been severely reduced by chemical recombination and plays 
a negligible role. To a plasma physicist this geometry is a classi¬ 
cal flute mode geometry and one might expect this equatorial iono¬ 
spheric geometry to be unstable to a variety of plasma instabilities. 
However, until recently, one of the basic difficulties was getting the 
unstable irregularities to the topside (i.e., above the F peak) when 
they were initiated on the bottomside. One must remember that the 
experimental evidence has exhibited both top and bottomside 
irregularities. 

At this point a general brief review will be given of equatorial 
Spread F (ESF) theories. First, we will discuss the linear theories. 
Dungey (1956) was the first to suggest that ESF was initiated on the 
bottomside by a Rayleigh-Taylor instability. Dagg (1957) suggested 
that the ESF phenomena was due to E to F region coupling, i.e., 
irregularities in the E region coupled up to the F region. In 1959, 
Martyn (1959) was the first to suggest that ESF was a manifestation of 
the E X B gradient drift instability. Calvert (1963) proposed that 
the downward motion of the neutral atmosphere at night was responsible 
for ESF. This mechanism is essentially equivalent to the E X B 






instability because of the relative motion between ions and neutrals 
in determining the instability. All of these previously invoked 
linear instability mechanisms could only explain the formation of 
bottomside irregularities. The collisional Rayleigh-Taylor instability 
with field line averaging was also proposed ( Balsley et al. , 1972; 
Haerendel , 1974) as a linear instability mechanism. By averaging 
(integrating) the density along the magnetic field, the total electron 
content profile becomes steeper on the bottomside and its peak is 
raised in altitude with respect to the local electron density peak. 

This would allow the linear mechanism to operate to slightly higher 
altitudes 100 km greater), but still would not explain the 
existence of irregularities above this "new peak". Hudson and 
Kennel (1975) pointed out the importance of the collisional drift 
mode in ESF in the wavelength regime 30m - 100m. This mode could be 
excited on both the top and bottomside but still would not explain the 
longer wavelengths. In their paper, finite larmor radius (FLR) 
corrections were also applied to the collisionless and collisional 
Rayleigh-Taylor instability. 

Several nonlinear theories have been invoked to explain the 
different ESF observations. For example, Hudson et al. (1973) sug¬ 
gested that the very smallest scale (*«.10m) irregularities (e.g., 
those seen by radar coherent backscatter) were due to a two step 
process. In this prescription a longer wavelength instability sets 
up the driving conditions for the shorter wavelengths to become 
unstable. This is similar in spirit to the successful two step 
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theory ( Sudan et al. , 1973) proposed for Type II equatorial E region 
electrojet irregularities. Haerendel (1974) suggested that the range 
of wavelengths (many kilometers down to centimeters) exhibited by ESF 
phenomena was due to a multi-step process. This scenario is as 
follows: (i) the collisional Rayleigh-Taylor (R-T) instability with 

horizontal wavevectors is driven by gravity and the background, zero 
order electron density gradient scale length on the bottomside; then 
(ii) the E X B gradient drift instability with vertical wavevectors 
arises due to the horizontal density, large amplitude variations set 
up by the collisional R-T instability; then (iii) the inertia 
(collisionless) dominated R-T instability arises; and finally (iv) 
kinetic drift waves grow upon these irregularities after they reach 

_2 

large amplitude. Chaturvedi and Kaw (1976) interpreted the k 
measured power spectrum of the ESF plasma density irregularities in 
terms of a two step theory. In this theory longer wavelength R-T 
modes couple to kinetic collisional drift waves in such a manner that 

_2 

the mode coupling results in the observed k spectrum. 

A major breakthrough was made by Scannapieco and Ossakow (1976) 
who performed a nonlinear numerical simulation of the collisional R-T 
instability for ESF geometry. The simulation results showed that the 
collisional R-T instability generated irregularities and bubbles 
(plasma density depletions) on the bottomside of the F region which 
subsequently rose beyond the F peak by nonlinear polarization induced 
E X B forces. This was the first theoretical result to explain how 
long wavelength irregularities could appear on both the bottomside and 







topside of the F region. The bubble phenomena was in accord with the 
recent observations ( Kelley et al. . 1976; McClure et al. , 1977, 

Woodman and La Hoz , 1976) of plasma density depletions. An analytical 
nonlinear mode-mode coupling theory for the coherent development of 
the collisional R-T instability was performed by Chaturvedi and 
Ossakow (1977). This theory suggested that vertical modes would be 
dominant and result in a k 2 power spectrum. Hudson (1978) extended 
the previous results to the collisionless R-T regime and reached 
similar conclusions. Analytical models for the rise of collisional 
and collisionless R-T ESF bubbles, in analogy with fluid bubbles, 
was presented by Ott (1978). At the same time, Ossakow and Chaturvedi 
(1978) presented analytical models for the rise of collisional R-T ESF 
bubbles within the context of the electrical analogy with barium 
clouds. 

Costa and Kelley (1978a,b) suggested that coherent steepened 
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structures and not turbulences would give a k power spectrum. More¬ 
over, these sharp gradients could cause small scale sizes (^ 20m) by 
collisionless low frequency (much less than the ion gyrofrequency, Q^) 
kinetic drift waves via a two step process. Their analysis was a 
linear one carried out on a nonlinear state, i.e., one achieves the 
steepened gradients by nonlinear processes and then one performs 
linear theory on this state. Kelley and Ott (1978) suggested that the 
ESF bubbles, in the collisionless R-T regime, generate a wake with 
vortices. They then applied two dimensional fluid turbulence theory 


to the model. This resulted in the development of turbulence at 





shorter and longer wavelengths than the bubble size. This in turn led 
to a prediction of k 1 for the power spectrum (which does not appear 
to be in agreement with existing experimental observations) in the 
range L g 1 < k < L D *, where L s the st * rri - n 9 (bubble) size and Lp 
is a dissipation length cutoff. In a continuation of the numerical 
simulation work, Ossakow et al. (1978) showed a more rapid ESF develop¬ 
ment and higher bubble rise velocities resulting from sharper bottom- 
side background electron density gradients and higher altitudes of the 
F peak. In Huba et al. (1978) very small scale (wavelengths 'v, lm and 
36cm) irregularities are reported. A two step process, utilizing high 
frequency kinetic drift cyclotron or lower hybrid drift 

instabilities, is invoked to explain them. Linear theory for these 
instabilities, in the ESF environment, was performed on the nonlinear 
state with encouraging results. 

The above introductory remarks glaringly point out that much work 
and significant progress in the theoretical area of equatorial Spread 
F has been accomplished in the past few years (indeed just look at 
the number of publications during 1978 alone). Notwithstanding the 
recent successes, much work still needs to be done. Indeed, the 
theoretical and numerical simulation efforts in ESF are continuing 
along a hot and heavy path. Section II of this paper presents out¬ 
lines of some of the theoretical efforts briefly mentioned in the 
preceding paragraphs. Given the length limitations, it would be 
exceedingly difficult to outline all of the theoretical works 


mentioned or even to present all of the details of a few works. 







Section II, hopefully, will wetten the reader's appetite to read the 


referenced works. Section III presents a summary concerning ESF 
theory. 


II. Theory 

In this section we present some representative theoretical and 
numerical simulation works with the appropriate references. 

a. General . The basic plasma fluid equations applicable to the 
equatorial Spread F ionosphere are as follows: 

3n 


3t 


• V • (n V ) = P - ' 


-T Vn - en (-V<(> + - 


m i n i ( 3t + -i * V) = " T < Vn < 


i i 2. i i in—i 


( 1 ) 


£ = ne (V. - V^) (5) 

In the above equations the subscript a denotes species (e is electron, 
i is ion), n is density, V is velocity, P is the production, v R is the 
chemical recombination rate, T is temperature, V is the gradient 











operator, the ambient magnetic field (taken to be uniform), 

e is the electronic charge, c is the speed of light, m is mass, £ is 
gravity, is the ion-neutral collision frequency, J is current, and 
the electrostatic approximation has been made where E - —V«J>. 

Equation (1) is the continuity equation, (2) and (3) are the electron 
and ion momentum equations, respectively, (4) is the divergence of 
the current and (5) is the current equation. What we have in mind is 
to apply the set of equations (1) - (5) to the two dimensions perpen¬ 


dicular to B at the geomagnetic equator, making various approxima¬ 


tions. 


Assuming a harmonic perturbation dependence of the form 













Making the electrostatic assumption 

E - - V* (12) 

and breaking the potential into an equilibrium and a perturbed 
quantity, 

♦ = 4> 0 * $ (13) 

eqns. (8) - (11) become 

f£ _ |-(V$ x z)V • n = - v R (n-n Q ) (14) 

B 

V • (v ±n n V^!) = ^( £ x z) • Vn (15) 

where eqns. (14) and (15) are taken to be two dimensional ( A 1^). 
Linearizing eqns. (14) and (15), taking a horizontal perturbation 
results, as in eqn. (7a) but in a more illustrative form, in the 
linear growth rate 



This clearly shows that only the bottomside of the F region where 
Vn Q is positive can be linearly unstable (and only if the first term 



Equations (14) and (15) were solved numerically using a 
vertical mesh spacing of Ay« 2 km and total y extent of 200km, and 
an east-west, horizontal mesh spacing Ax - 200m and a total horizontal 
extent of 8km. Realistic profiles of v^ n and v R as a function of 
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altitude were utilized. The system was initialized with a 

perturbation of a few percent in the horizontal (x), east-west 

direction with a wavelength , v3km and the evolution in time of (14) 

and (15) was followed for different background electron density 

profiles. Figure 3 shows the results for a background electron 

density, n Q , profile with an F peak at 354km and a minimum bottomside 

background electron density gradient scale length, L , ''10km. In this 

case at t * 4000 sec a bubble (plasma density depletion) is clearly 

forming and beginning to rise in the central portion of the mesh 

(note n i n + n,). The isodensity contours are such that the maximum 
o x 

absolute value of the enhancement or depletion is in the center and 
the contours decrease (in absolute percentage) as one goes toward the 
outer contours. At t = 4000 sec the maximum depletion within the 
bubble is 54% and the maximum enhancement over the mesh is 84%. At 
t = 8000 sec we notice that the bubble has reached the altitude of 
the F peak, with the innermost contour of the rising bubble repre- 

4 

senting a 41% depletion. At t = 10 sec the main bubble is clearly 
through the F peak with an innermost depletion contour of 41%. How¬ 
ever, in the ionosphere below the bubble near x = 0 there is a 71% 
depletion contour, similarly in the wings near |x| * 4km. The inner¬ 
most enhancement contour, at this time, represents a 236% enhancement 
with a maximum inside this contour of 294%. Note that the top of the 
main bubble is at an altitude of 375km while the bottom trail of the 

4 

bubble is at an altitude of 270km. Between t = 8000 ard 10 sec the 
bubble has risen ''-24km which represents a rise velocity '''12m/sec. 
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Also note that the bubble is ''dkm wide. Figure 4 depicts contours 

4 

of constant induced potential ^ at t - 10 sec. This shows that 
the more isolated part of the high altitude bubble depicted in Fig. 3 
is acted on by an induced electric field which points from west to 
east and is dipolar in nature. This causes the bubble to rise with a 
(-c/B q )V ^ x z velocity. However, the lower portion of the mesh is 
acted on by an induced electric field v»hich points from east to west. 
This field is much weaker than the induced field acting on the 
isolated portion of the central bubble. The lower altitude electric 
field causes the enhancements and depletions to move downward. Thus, 
the lower altitude portion of the central bubble becomes captured by 
the enhancements 

Figure 5 presents for comparison a case in which the background 
electron density profile (shape) was kept the same, but the entire 
profile moved up in altitude so that the F peak was at 434km. All 
other parameters are the same as in Figure 3, except and v R are 
taken for the altitude range 332km to 532km (those used in this 
simulation). One can immediately note the more rapid time evolution 
of the Spread F process with respect to that presented in Fig. 3. 

At t = 700 sec, a rising bubble with an innermost contour of 79% and 
a maximum depletion inside this contour of 84% were noted. At 
t = 1000 sec the bubble has reached the peak and at this time the 
innermost depletion contour is 85%. The long trail associated with 
the high altitude bubble has a 100km extension to lower altitudes. 

At t * 1400 sec the top of the main bubble is at an altitude 'vSOOkm 










and has a long trail connecting to an altitude of 357km. There is a 
maximum 70% depletion within the innermost contour of the high 
altitude bubble. Between t = 10 3 and 1400 sec the top part of the 
bubble rose ' v '65km and this represents a rise velocity ^leOm/sec. 
Potential contour results for this simulation show similar patterns to 
those exhibited in Fig. 4. Naturally, the induced electric fields 
causing the bubble to rise in the present case are stronger. This 
spread in bubble rise velocities has been observed by AE satellite 
data ( McClure et al. , 1977). 

Other numerical simulations in this series have been performed 
and the paper by Ossakow et al. (1978) should be consulted for more 
details. The basic conclusions reached from these simulations are as 
follows: (i) the collisional R-T instability causes linear growth on 

the bottomside of the equatorial Spread F region; (ii) plasma density 
depletions (bubbles) steepen on their top and nonlinearly rise to the 
topside ionosphere, beyond the F peak, by polarization (induced) 

E X B forces; and (iii) high altitude of the F peak, small bottomside 
background electron density gradient scale lengths, and large percent¬ 
age depletions yield large vertical bubble rise velocities, with the 
first two conditions favoring collisional R-T linear growth (instabil¬ 
ity). In addition, large spatial bubbles with similar rise velocities 
to those presented here, but with almost 100% depletions, have been 
produced by numerical simulations ( Zalesak et al. , 1978). In these 
cases the horizontal mesh covered 200km in extent (intermesh spacing 
Ax = 5km) and a long wavelength ('WSkm) initial perturbation was 
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used. Bubbles with horizontal dimensions '\'50km resulted. The NRL 


group has also added neutral wind effects to these simulations and 
found that an eastward neutral wind results in a westward motion of 
the bubbles in addition to its rise. This is also in agreement with 
observations of bubble motion (see McClure et al. , 1977). 


c. Analytical 2D Coherent Mode Coupling Results. A two 



dimensional nonlinear quasi-final state of the collisional R-T 


instability was investigated by Chaturvedi and Ossakow (1977) using 


analytical means and considering coherent mode coupling as the 


saturation mechanism, 


etc. This yields the following coupled nonlinear 


equations 


the first 


term on the RHS represents damping, and the second term 


the nonlinear term. Comparing the nonlinear term in (16) with the 


last term in (17) results in 


Therefore, eqn. (17) is treated linearly and the nonlinearity is 


retained in (16). A perturbation of the form 
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For typical values of the parameters, A. 


(i.e., the power spectrum would = k 


coherent nonlinear evolution where the linearly damped mode A. 


generated nonlinearly by the linearly unstable mode A 


contains 2k ). More 


contains k and A. 


detailed information can be found in Chaturvedi and Ossakow (1977) 


d. Analytical Models for ESF Bubble! 


.11 discuss 


collisional R-T bubbles ( Ossakow and Chaturvedi , 1978) in which the 


electrical analogy with plasma density enhancement (e.g., barium 


clouds) was utilized and results obtained for general 2D (Jfi^) 


bubble shapes. Here eqns. (14) and (15) 


itilized with the 


further simplifying assumptions of neglecting recombination chemistry 


•neutral 


0) and the explicit altitude dependence of the 


The following set of equations result, 


collision frequency, 


where an ambient horizontal electric field E has been included to 


show more generality. Equation (24) can be thought of, in general 


as a potential equation for a dielectric immersed in the applied 


The plasma density depletion (bubble) is 


electric field, 


analogous to the case of a cavity immersed in a dielectric with a 









uniform electric field (E*). In general (23) and (24) do not admit of 
a two dimensional (jB q ), steady state solution. However, for the case 
of a constant density inside the depletion and a constant density out¬ 
side the depletion, two dimensional steady state solutions can be 
obtained. 

For a constant density depletion with an elliptical shape one has 
n(x,y) * n Q - rip^y) 


"d 


■ "o f 1 - ir F(x '*>] 

fin H [l (|) 2 - ( r) 2 J 

11. x > 0 

H (x) = ] 

to, x < 0 


(26) 


(27) 


where H is the Heaviside function (note the geometry is taken such 
that x is vertical and y horizontal). For this elliptical shape, 
the solution to (24), neglecting E , is 


it 


3y 


°i V / 3 5n/n o \ 

V in e [b* a(l - 6n/n o )J 


(28) 


Using -(c/B q ) 7$ x z this further results in the nonlinear vertical 
bubble rise velocity, V g , given by 

g / a 6n / n o \ 

V B*v in ( b + a( l-6n/n o )J < 29 > 

Limiting cases of (29) for sheet, cylindrical and slab bubble 
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geometries are given by 



(6n/n Q ) 

(1 - 

6n/n o ) _1 , 

b « a 

(30a) 

(6n/n Q ) 

(2 - 

6n/n o ) -1 , 

b = a 

(30b) 


0 

, 

b << a 

(30c) 


For typical values of v in as a function of altitude, Fig. 6 exhibits 



rise velocities given by (30b) for various percentage depletions. 
Table 1 shows the rise velocities, in units of g/v^ n , for various 
shapes. (Note that the linear case comes from linearizing (23) and 


(24)). All of the geometry results can be expressed in a concise 
formula 



where f(Sn/n Q ) is an increasing function of the percentage 
depletion 6n/n Q . Basically, the results predict that high altitudes 
and/or large percentage depletions yield high vertical rise velocities 
for the bubbles (in agreement with experimental observations). 

Collisional and collisionless (inertia dominated) two dimensional 
cylindrical R-T bubbles have been studied analytically by Ott (1978). 
This study is based on the analogy with fluid dynamic flows and brings 
forth some of the work done on two-dimensional fluids. This study 
begins with the basic equations (1) - (3), considers two dimensions 
(ABq) , sets n fi as n^ n, and makes the assumption that 


To lowest order, using (31), one obtains from (3) a lqwest order ion 
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velocity (with z 


(Note that inks units are being used here to coincide with the units 


used by Ott (1978)). Quasi-neutrality, i.e., 7 • J = 0, with the 


assumption of two dimensionality implies that J can be specified in 


terms of a single scalar potential function, <(», such that 


Using eqn. (2) for the electron velocity, and a next order ion 


velocity equation obtained by putting (32) into (3), the following 


ion velocity equation is obtained 


with P ■ P + p. + B i|». Using (32) and the assumption af either 


isothermal or adiabatic ions one has 


in the ion continuity equation (1) with v 


Using V, 






Equations (34) - (36) form a complete set of equations 
sufficient to determine the unknown quantities , n, and p. In 

the limit of v^ n -*■ o, eqns. (34) - (36) are identical to those of 
an ideal incompressible fluid. At this point the philosophy taken 
is that there is much to be learned concerning bubbles in the ESF 
ionosphere from the extensive studies of bubbles in fluids. One 
then uses a stream function £ such that 

V. (o) s i x V ? (37) 

After a series of manipulations one can obtain an equation for the 
stream function. In the collisional and collisionless case one 
finds that for certain values of the* bubble rise velocity £ will 
satisfy the equation, for a cylindrical shape (see Ott , 1978 for 
more details). 

The results of this study by Ott (1978) predicts the bubbles 
to be cylindrical (circular cap at top) in two dimensions. The 
bubble rise velocity in the collision-dominated regime, for a 100% 
depletion, is given by 

V = -— (38) 

B v. 
in 

which is altitude dependent. This result is the same as that pre¬ 
dicted by (30b) for 6n/n Q = 1, i.e., a 100% depletion. In the 
inertia (collisionless) dominated regime, for a 100% depletion, the 
bubble rise velocity is given by 

V B * «j (Rg) ** (39) 
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where R is the radius of curvature at the top of the bubble so that 
(39) is size dependent. 

e. Analytical Work on the Very Small Scale (^10m) 
Irregularities . The work in this section represents essentially 
multilinear calculations using kinetic theory resulting in plasma 
kinetic drift wave modes. It is multilinear because it depends on 
say a two-step process whereby linear theory is performed on the non¬ 
linear state. The driving density gradients, in these calculations, 
are thought to arise through a primary instability, driven by the zer 
order background ionospheric equatorial F region electron density 
gradient, achieving a large amplitude state. The zero order iono¬ 
spheric electron density gradient is of larger scale size than the 
primary instability electron density gradient scale size, which would 
be of the order of the instability wavelength. Because the calcula¬ 
tions are kinetic, they employ particle distribution functions. 
Kinetic drift waves have been investigated for laboratory plasma 
fusion conditions for over twenty years, so a well developed formalis 
exists. 

Before proceeding with specific calculations let us present some 
general concepts regarding kinetic drift waves which will be useful 
for both types of calculations presented in this section. The basic 


geometry is such that 










(similar to the zero order background equatorial ionospheric 
geometry) where here n Q (x) arises due to the primary instability. The 
orthogonal coordinate system is completed by the ion and electron 
diamagnetic velocities being along the y axis, 


^d = (V di - V de ) y 


V = (V. /2fi.) (3ln n /3x) 


V = - (V /2£2 ) (31n n /3x) 


where g = (2^ e /nK g ) and the larmor radius is defined by 

r a = v a / fi a - The linear analysis is then performed with this in mind. 

First the low frequency (u> « ft^) collisionless drift wave 
calculations of Costa and Kelley (1978a,b) , for k^ 'v* 1 (r^ is the ion 
gyroradius), will be presented. These calculations were meant to pro¬ 
vide a basis within which to try and account for the 3 meter radar 
backscatter observations of Jicamarca (see Woodman and La Hoz , 1976). 


A linear kinetic dispersion relation is derived using perturbations of 

the form exp i k • x - wt , where k = k.. z + k y and u> = u + iy. 

L J ]| j. r 

Furthermore, assumptions are made such that ui « ft., V. <<|u>/k,, |« V 

ii II e 

and u > v^ n# v^ g , v^, which are the ion-neutra^ ion-electron, and 
ion-ion collision frequencies, respectively. This analysis depends on 
having a k ^ , a component of the wavevector along the ambient geo¬ 

magnetic field. Figure 7 shows some growth rate results from these 
calculations. The growth rate is in units of the ion thermal velocity 
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divided by the electron density gradient scale length 
(L = n Q (dn Q /dx) *). Figure 8 depicts some measured inverse electron 
density gradient scale lengths during ESF and maximum growth rates as 
a function of these inverse scale lengths. The basic results of these 
calculations show maximum growth for k r^ « 1.5, i.e., \ » 20m (for 

typical ESF parameters) with growth rates ^ 1 sec 1 . For more 
detailed analysis concerning these low frequency drift waves applied 
to ESF see Costa and Kelley (1978a,b). 

Now we present the high frequency drift wave analysis of Huba 
et al. (1978). In this reference radar backscatter observed 
irregularities with wavelengths of 1 meter and 36 cm at Kwajalein dur¬ 
ing ESF conditions are shown. In an effort to explain these very short 
wavelength irregularities high frequency (o> ^ (h) drift waves were 
analyzed for ESF conditions. These waves are the so-called drift 
cyclotron (DC) or lower hybrid drift (LHD) instabilities with maximum 
growth rates for k r ^ 1 (and k. = 0) . No k.. is required for 

A 6 II ll 

these instabilities. The parameter determining which instability 
operates in a collisional plasma is 

C f = (v ii /fl i ) (k x r i )2 (42) 

Utilizing the linear dispersion relation for high frequency drift 
waves for C^ << 1, instability results for 

L/r. < (1/2*) (m./m J* 5 (43) 

l i e 

where * is the harmonic number (o^ « *J2^). For the 0 + ESF ionospheric 
plasma this requires the electron density gradient scale length 







L < 340m, which is satisfied (see Fig. 8). Growth rates for 
these instabilities are given by 


s (m /m.) l «. 


(44) 


However, the condition C„ << 1 implies that (k r.) 2 n << 2 x 10 7 and 
f x 1 

for k r g -v l this means n « 10 3 cm 3 , which is quite restrictive. 

Longer wavelengths, i.e., smaller values of k r would raise the 

x e 

density restriction somewhat, but still be restrictive for ESF 
conditions. 

For C f ^ 1, the lower hybrid drift instability is operative and 
there is no threshhold condition on L. Basically, the collisions 
which increase C f destroy the ion gyroresonances needed for the DC 
instability to operate. The real and imaginary part of the fre¬ 
quency for the LHD instability are given by 


(?) 


(45) 


In the collisionless limit it should be noted that the DC instability 
transforms into the LHD instability for high enough ion diamagnetic 
velocities such that 
L/r. 6, (m./m )^ 


which for ESF conditions implies that L ^ 30m. Figures 9 and 10 show 
some typical results from the analysis of the high frequency drift 


wave linear dispersion relation. 








The results of this analysis predicts that the lower hybrid 
drift instability is dominant for most typical ESF ionospheric para¬ 
meters. Also, maximum growth of the instability occurs for 

k r -v, 1 (X ^21 cm), although good growth rates can occur for 
u e j. 

X -v. 1 m. Finally, from this instability, large growth rates 
(y £h) resulting in growth times, t = y *» less than a second can 
occur. For more details of this work see Huba et al. (1978). 

III. Summary 

Although much progress has been made in the theoretical efforts 
directed toward the equatorial Spread F ionosphere, especially in the 
past three years, more has to be done. Also the burden cannot be 
placed on the theoretician alone. Correlative measurements have to 
be made prior to and during ESF conditions. It is not sufficient to 
make a single measurement with one instrument and then expect a 
complete theoretical description of ESF. One needs to know the state 
of the ionosphere with respect to driving parameters such as back¬ 
ground electron density profiles (bottomside electron density 
gradient scale lengths and height of the F peak), d.c. electric 
fields, neutral winds, and ionic mass composition in order to build a 
predictive model. In addition, in order to compare results from the 
predictive model, measurements of the in situ fluctuating component 
of the electric field and plasma density have to be made, as well as 
radar backscatter measurements of the very small scale irregularities 
(V10m) and ground measurements of satellite signal propagation 
amplitude and phase scintillations. 
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Soph? of the remaining theoretical problems could be listed as 
follows. (1) What are the effects of changing initial conditions, 
including ion inertia, and including neutral winds in the numerical 
simulations? (2) An analytical description of many bubbles rising. 

(3) How do bubbles decay and what role does diffusion, etc. play? 

(4) How does a turbulent development occur? (5) What are the effects 
of other regions of the ionosphere (e.g., E region) on ESF? (6) A 
more complete study of collisional effects on drift waves is needed. 
What determines when the small scale irregularities should occur (a 
more quantitative description)? (7) A determination of the nonlinear 
saturation of small scale irregularities (instabilities) is needed. 

(8) What are the effects of k ? indeed points (6) and (7) are tied 
to the more general question of what is the relation between the very 

small scale (6,10m) irregularities (which radar backscatter observes) 
and the longer wavelength fluid type (e.g., R-T) irregularities (which 
are primarily responsible for Spread F seen on ionograms)? I am sure 
that some more questions and points could be raised. However, the 
above list should keep the theoreticians busy for a reasonable time. 
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Fig. 3 - Contour plots of constant ni/n 0 for the simulation with an F 
peak at 354km at t * 2000, 4000, 8000, and 10,000 sec. The smell dashed 
contours with a plus sign inside and the solid contours with a minus sign 
inside indicate enhancement and depletions over the ambient electron 
number density. The large dashed curve depicts the ambient electron 
number density (values on upper horizontal axis), n Q , as a function of 
altitude. The vertical y axis represents altitude, the lower horizontal 
x axis is east-west range, and the ambient magnetic field is along the 
z axis, out of the figure. Taken from Ossakow et al. (1978) 














Fig. 4 - Contours of constant induced (polarization) potential,, 
over the mesh (see Fig. 3) at t ■ 10,000 sec. Plus and minus denote 
positive and negative values, with values decreasing in magnitude as 
one goes from the innermost to the outermost contours. The large 
dashed curve ia n c . Taken from Ossakow et al. (1978) 










Fig. 5 - Contour plots of constant nj/n*, for the F peak at 
434km simulation at t - 300, 700, 1000, and 1400 sec. All 
other nomenclature is same as Fig. 3. Taken from Ossakow 

























PERPENDICULAR PARAMETER, b = k JL /3 i /2 

Fig. 7 - Growth rates as a function of b for the low frequency collisionless drift waves 
(Pi is the ion gyroradius). The curve Y << u r is obtained from a small growth rate appro: 
imatlon expression. Taken from Costa and Kelley (1978b) 










JAVELIN 8:63 NOVEMBER 18, 1973 NATAL, BRAZIL 



Fig. 8 - Observed Inverse gradient scale lengths (1/n) (dn/dx) as a function of altitude 
during the Natal rocket flight. Also indicated by the straight line is the maximum growth 
rate cf low frequency drift waves as a function of (1/n) (dn/dx). Note that p. is the ion 
gyroradius. Taken from Costa and Kelley (1978a) 1 
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Fig. 10 - Maximum growth rate as a function of ion diamagnetic 
drift velocity for the DC instability (solid line) and the LHD 
instability (dashed line). The other parameters are the same 
as in Fig. 9. Note that the DC instability goes into the LHD 
instbility for Vjj/Vi > 0.11. Taken from Huba et al. (1978) 
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